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We calculate the charge form factor and the longitudinal structure function for 16 and 
compare with the available experimental data, up to a momentum transfer of 4 fm _1 . The 
ground state correlations are generated using the coupled cluster [exp(S)] method, together 
with the realistic vl8 NN interaction and the Urbana IX three-nucleon interaction. Center- 
of-mass corrections are dealt with by adding a center-of-mass Hamiltonian to the usual 
internal Hamiltonian, and by means of a many-body expansion for the computation of the 
observables measured in the center-of-mass system. 
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One of the fundamental problems in nuclear physics is 
related to developing a complete understanding of how 
nuclear structure arises as a result of the underlying in- 
teraction between nucleons. This in turn should help 
us develop a complete understanding of the electromag- 
netic structure of the nucleus, as revealed by the wealth 
of high-quality data that electron scattering experiments 
have provided for the past 30 years. The interplay of 
nuclear correlations, meson-exchange current or charge 
densities, relativistic effects in nuclei, the importance of 
three- or more many-body interactions in relation with 
the dominant two-body interaction in nuclei awaits to be- 
ing assessed in greater detail. Unfortunately, solutions of 
the many-body Schrodinger equation with realistic inter- 
actions have been proven very difficult to obtain. Only 
in recent years has progress been made and first results 
of microscopic calculations relating to ground-state and 
low-excited states for nuclei with A < 7 been reported M . 
These calculations have been obtained using the Green's 
function Monte Carlo method, but this approach, just 
like the Fadeev H or the Correlated Hyperspherical Har- 
monics H methods successfully used for the A = 3,4- 
body system, suffers limitations in the number of nucle- 
ons they can treat. To date, only the Variational Monte 
Carlo method |4| has enjoyed success in solving the many- 
body problem for medium nuclei, but those results still 
show a room for improvement. 

We are using the exp(S) coupled-cluster expansion to 
calculate the ground state of 16 0. Our approach is very 
similar to the standard approach, first developed by the 
Bochum group [g , and has been outlined recently in [|S| . 



The idea behind the coupled-cluster expansion formal- 
ism relies on the ability of expanding the model nuclear 
wave function in the many-body Hilbert space in terms 
of two Abelian subalgebras of multiconfigurational cre- 
ation and their Hermitian-adjoint destruction operators. 
The expansion coefficients carry then the interpretation 
of nuclear correlations. The fact that we make no artifi- 
cial separation between "short-range" and "long-range" 
correlations is one particular strength of this many-body 
method. 

The derivation of the explicit equations is quite te- 
dious, but requires only standard techniques. For a 
closed-shell nuclear system, the total Hamiltonian is 
given as 
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The Hamiltonian includes a nonrelativistic one-body ki- 
netic energy, a two-nucleon potential, and a supplemen- 
tal three-nucleon potential. We have chosen the Argonne 
wl8 potential |6|] as the most realistic nucleon-nucleon in- 
teraction available today. The Argonne vl8 model pro- 
vides an accurate fit for both pp and nn scattering data 
up to 350 MeV with a x 2 /datum near one. The introduc- 
tion of charge-independence breaking in the strong force 
is the key element of obtaining this high performance. 
However, the two-body part of this interaction results in 
over-binding and too large a saturation density in nuclear 
matter. Therefore, the NN potential is supplemented 
by a three-nucleon interaction (part of the Urbana fam- 
ily n), which includes a long-range two-pion exchange 
and a short-range phenomenological component. The 
Urbana-IX potential is adjusted to reproduce the bind- 
ing energy of 3 H and give reasonable saturation density 
in nuclear matter when used with Argonne i>18 [fij. 



We are searching for the correlated ground state of the 
Hamiltonian H, which we denote by |0). The ansatz for 
the many-body wave function 0) is defined as the result 
of the cluster correlation operator, S^, acting on the ref- 
erence state of the many-body system, the uncorrelated 
ground state |0): 
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For a number-conserving Fermi system, the standard 
choice for |0) is the single-particle shell-model (Slater 
determinant) state formed from an antisymmetrized 
product of single-particle wave functions. The clus- 
ter correlation operator is defined in terms of its ph- 
creation operators expansion (Oj = 1, Oj = at aj^, 

°2 = 4i42 a fc2 a /n) as 
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The problem of solving for the many-body wave function 
|0) and the ground-state energy, E, is now reduced to 
the problem of solving for the amplitudes S n . This im- 
plies solving a set of non-linear equations, which may be 
obtained using a variational principle. We construct a 
variation S\0) orthogonal to the correlated ground state 
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and require that the Hamiltonian between the ground 
state and such a variation vanishes. As a result, we ob- 
tain an equation for the ground-state energy eigenvalue 
E in terms of the cluster correlation coefficients, {S n }, 
and a set of formally exact coupled nonlinear equations 
for these coefficients: 
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Then, the computation breaks down into two steps |g|: 
In the first step, the G-matrix interaction is calculated 
inside the nucleus including all the corrections. This re- 
sults in amplitudes for the 2p2h correlations, which are 
implicitly corrected for the presence of 3p3/i and ApAh 
correlations. In the second step the mean field is cal- 
culated from these correlations and the single-particle 
Hamiltonian is solved to give mean-field eigenfunctions 
and single-particle energies. These two steps are iter- 
ated until a stable solution is obtained. Calculations are 
carried out entirely in configuration space where a 50fkv 
space is used. The general approach, when the Hamil- 
tonian includes only up to two-body operators, has been 
presented in MM. The results we report here have been 
obtained by taking into account the three-nucleon inter- 
action via a density-dependent approach. The details of 



Once the correlated ground state |0) is obtained, we 
can calculate the expectation value of any arbitrary op- 
erator A as 

o= {0|e s Ae- s S t |0) , 

where S^ is also defined by its decomposition in terms of 
p/i-creation operators 

S = Z^ ~\ Sn °™ • 



The amplitudes S n are obtained in terms of the S n am- 
plitudes in an iterative fashion. 

Note that the correlated ground state |0) is not transla- 
tionally invariant since it depends on the 3A coordinates 
of the nucleons in the laboratory frame. Therefore, in 
practice one has to take special care of correcting for the 
effects of the center-of-mass motion. This is done in sev- 
eral steps: First, the Hamiltonian (jjj) is replaced by the 
internal Hamiltonian 
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which is now entirely written in the center-of-mass frame 
by removing the center-of-mass kinetic energy, Tqm — 
P^ M /(2mA), with m the nucleon mass. Both the two- 
and three-nucleon interactions are given in terms of the 
relative distances between nucleons, so in this respect no 
corrections are needed. Secondly, a many-body expan- 
sion has been devised pij| in order to carry out the nec- 
essary corrections required by the calculation of observ- 
ablcs, which are measured experimentally in the center- 
of-mass frame. This procedure is based on the assump- 
tion that we can neglect the correlations between the 
center-of-mass and relative coordinates degrees of free- 
dom, and a factorization of the correlated ground state 
|0) into components which depend only on the center-of- 
mass and the relative coordinates, respectively, is possi- 
ble. We also assume that, indeed, the correlated ground 
state |0) provides a good description of the internal struc- 
ture of the nucleus. Finally, in order to ensure such a 
separation, a supplemental center-of-mass Hamiltonian 
is added 
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which has the role of constraining the center-of-mass 
component of the ground-state wave function []12j. We 
choose the values of the parameters (3cm and fl 2 such 
that they correspond to a value domain for which the 
binding energy 

E = {Hint = Hint + H CM ) - (H C m) 

is relatively insensitive to the choice of the Pcm and fi 2 



this approach will be presented elsewhere |10| 



values 1 1 1 . When leaving out the center-of-mass Hamil- 
tonian, the calculated binding energy of 16 is equal 



to 7.54 Mev/nucleon, which is thought as a reasonable 
value, given the uncertainties related to the three-nucleon 
interaction. 

Figure |l] shows the theoretical result for the charge 
form factor in 16 0. In the one-body Born-approximation 
picture, the charge form factor is given as 

F L (q) = (6\Y,Mq 2 )e^\o), 

k 

with fk(q 2 ) the nucleon form factor, which takes into ac- 
count the finite size of the nucleon k |l3| . We also take 
into account the model-independent part of the "Helsinki 
meson-exchange model" |14| , namely the contributions 
from the %- and p-exchange "seagull" diagrams, with the 
pion- and p-meson propagators replaced by the Fourier 
transforms of the isospin-dependent spin-spin and ten- 
sor components of the ^18 NN interaction. This sub- 
stitution is required in order for the exchange current 
operator to satisfy the continuity equation together with 
the interaction model. The contributions of the n- and 
p-exchange charge density give a measurable correction 
only for q > 2 fin -1 . 

In order to generate the form factor depicted in Fig. [jL 
we have first used the procedures of pl[ , keeping the 
contributions that can be written in terms of the one- and 
two-body densities. Then we take the Fourier transform 
in order to produce the theoretical charge density. Using 
this theoretical charge density, we generate the charge 
form factor in a Distorted Wave Born Approximation 
picture (l5[| , in order to take into account the distortions 
due to the interaction of the electron with the Coulomb 
field. This last step results in smoothing out the sharp 
diffraction minima usually seen in the calculated charge 
form factor flM- The agreement with the experiment 
is reasonably good over the whole range of q spanned by 
the available experimental data ]lq| . 

A second electron scattering observable that we would 
like to compare with is the longitudinal structure func- 
tion SL(q), sometimes called the Coulomb sum rule, 
which is sensitive to the short-range correlations induced 
by the repulsive core of the NN interaction |17],[l|] . The 
Coulomb sum rule, SL(q), represents the total integrated 
strength of the longitudinal response function measured 
in inclusive electron scattering. In the nonrelativistic 
limit [ |L8| , we have 

S L (q) EE l + pLL ( q )-±\(0\p( q )\d)\^ 

where p(q) is the nuclear charge operator 



P(l) = \ E^'tlH.i), 



and pLh{q) is the longitudinal-longitudinal distribution 
function 



PLL(q) = d?i df 2 Jo(q\n -r 2 \) p (p,p) (ri,r 2 ) . 
Here p^'P' (r 1,7*2) is the proton-proton two-body density 

P^ p Hn,r 2 ) 

= 7 E< 8 I *(* - n)S(r 2 - r 3 )(1 + r z>i )(l + r z>j ) |6) , 



normalized as 
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In light nuclei reasonable agreement between the- 
ory and experiment is obtained for the Coulomb sum 
rule ]19| . In heavier nuclei, however, the experimental 
situation is a lot more controversial, since both a cer- 
tain lack of strength has been reported and because of 
the inherent difficulty of separating the longitudinal and 
transverse contributions in the cross section due to the 
distortion effects of the electron waves in the nuclear 
Coulomb field. Figure @ shows the calculated Coulomb 
sum in 16 0. Since no experimental data are available for 
16 O, we compare the results of the present calculation 
with the 12 C experimental data from pO] with an esti- 
mate |21[1 for contributions from large to. The large error 
bars on the experimental data are largely due to system- 
atic uncertainties associated with tail contribution j22[ . 
Preliminary theoretical results for 12 C are also shown and 
appear to follow closely the results for theoretical curve 
for 16 0. 

This calculation represents the most detailed calcula- 
tion available today, using the coupled cluster expansion, 
for a nuclear system with A > 8. This also represents a 
contribution to the on-going effort of carrying out mi- 
croscopic calculations that directly produce nuclear shell 
structure from realistic nuclear interactions. Similar cal- 
culations for other closed-shell nuclei in the p- and sd- 
shell are currently under way. 
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FIG. 1. Charge form factor for 0-16, obtained in the 
DWBA picture (v 18 + UIX + ME) 




FIG. 2. Coulomb sum for C-12 and 0-16, compared with 
"experimental" C-12 data which include theoretically deter- 
mined high-energy tail corrections 



